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Magnetically-trapped atoms in Bose-Einstein condensates are spin polarized. Since the magnetic 
field is inhomogeneous, the atoms aquire Berry phases of the Aharonov-Bohm type during adiabatic 
motion. In the presence of an eletric field there is an additional Aharonov-Casher effect. Taking into 
account the limitations on the strength of the electric fields due to the polarizability of the atoms, 
we investigate the extent to which these effects can be used to induce rotation in a Bose-Einstein 
condensate. 



PACS numbers: 72.15.R, 03.75.F, 05.30. J, 32.60 
I. INTRODUCTION 

The ability to control and manipulate clouds of 
trapped cold atoms has resulted in a series of interest- 
ing experiments on Bose-Einstein condensates. A very 
striking example is the creation and observation of quan- 
tized vortex states, which are expected to occur in rotat- 
ing atomic condensates , much in the same way as in 
liquid helium II. 

In this paper we study the quantum mechanical ef- 
fects related to the adiabatic motion of the spin-polarized 
atoms in strong external electromagnetic fields. As will 
be explained in section II below, "adiabatic" here refers 
to the dynamics of the spin; spin flip transitions are as- 
sumed to be suppressed by a large Zeeman gap. 

One of our main conclusions is that the magnitude of 
the electric field E required to induce a vortex state is 
very high and approximately given by the simple relation 



eER. 



0.5 MeV, 



(1) 



where R is the size of the condensate, and m e is the elec- 
tron mass. Since E is limited both by direct experimental 
difficulties, and by polarization effects of the atoms (that 
essentially change the trapping potential), this condition 
cannot be satisfied in present experiments. It is, how- 
ever, not excluded that future experiments with much 
larger condensates could be in a regime where these ef- 
fects would become important. Apart from this possible 
application, we hope that the theoretical methods used 
in this paper can be of some general interest. 

In the presence of magnetic and electric fields, the 
atoms aquire quantum phases during adiabatic motion, 
and these effects can, as originally shown by Berry, be 
represented by effective gauge potentials. The most fa- 
mous - and here also the most obvious - example is 
the Berry phase due to the adiabatic rotation of the 
spin along with an external magnetic field B = SB. 



For a closed path C the phase 7b aquired is given by 
h"fs = —hm s £l where hm s is the projection of the total 
spin along the z axis and £1 is the solid angle swept out by 
B(r) along the curve C. As shown by Berry, this phase 
can also be expressed as a line integral of a vector poten- 
tial, slb, or, using Stokes theorem, as a surface integral 
of the corresponding effective magnetic field, b^, 



h-lB = f dr ■ a_g = dS bB 



(2) 



In the following, we will refer to &b and as the Berry 
potential and Berry field, respectively, and to avoid con- 
fusion in the notation, we shall denote real external elec- 
tromagnetic fields and vector potentials with capital let- 
ters, and effective fields and potentials with small letters. 

Below we shall incorporate the Berry phase in the de- 
scription of the Bose-Einstein condensate by coupling a 
complex scalar condensate wave function to the vector 
potential ag. By noting that is a monopole field in 
the parameter space spanned by B(r), it is easy to find 
the corresponding expression in r space, (see, e.g., || or 

1) 

(b B )i = (V x a B )i = ~hm s e ijk (djB x d k B) ■ B, (3) 

Because of the monopole character of the field the corre- 
sponding vector potential has a string like singularity or, 
alternatively, is given by a nontrival fibre bundle. This 
complication will be of no relevance for this paper since 
we shall not attempt to explicitly solve the field equations 
in the presence of &b- 

In a geometry which is such that the particles move 
in a region where the effective magnetic field vanishes, 
there can still be a purely topological effect of the 
Aharonov-Bohm type, but more typically hs is nonva- 
nishing throughout space and the particle experiences a 
Lorentz force (which is much smaller than the "Stern- 
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Gerlach" force V(/J-B) responsible for the magnetic trap- 
ping of the atoms). 

When electric fields are present, new effects occur. 
Qualitatively, this is because an atom moving with ve- 
locity v in an electric field E, will in its rest frame see a 
magnetic field B' given by 



B' = B-^vxE + C(v 2 /c 2 ), 



(4) 



and get a Zeeman shift oc jl ■ B' where pi is the magnetic 
moment of the spin polarized atom. This effect can be 
incorporated in the Lagrangian describing the system by 
an effective vector potential 



&AC = fix E/c 2 



(5) 



which is the origin of the Aharonov-Casher (AC) effect 
0. Just as in the case of the Berry field strength dis- 
cussed above, the Aharonov-Casher field, hAc, is nonva- 
nishing for a general configuration of the external electric 
and magnetic fields. In particular, this can be the case 
for a constant B field, as long as E is space dependent, 
or vice versa. The latter case is the most interesting for 
magnetically-trapped atoms. 

Effective magnetic fields in Bose-Einstein condensates 
can induce rotations, just as an ordinary B field in a type 
II superconductor. However, for realistic geometries the 
Berry potential can only give a phase < 2irm s , and 
since a single vortex corresponds to a 2ir rotation of the 
phase of the condensate wavefunction, this can at most 
induce a couple of vortices in condensates of alkali atoms. 
Since the strength of the AC potential is oc E, it can in 
principle give large phases, but since it is a relativistic 
effect, this would require very strong electric fields. An 
intersting feature of the AC potential is that the cor- 
responding effective magnetic field, hAc has a spatial 
distribution very different from that of the effective AB 
magnetic field. 

The physical origin of Berry and Aharonov-Casher 
phases are well understood, and the aim of this paper 
is to investigate whether they are of importance in real- 
istic configurations of atomic Bose-Einstein condensates. 
However, to make the presentation more self-contained, 
Sec. II contains a short summary of how the effective po- 
tentials &b and slac can be derived, and also how they 
can be incorporated into the Gross-Pitaevskii equation 
satisfied by the condensate wave function. In Sec. Ill 
we estimate the magnitude of the induced phases, and 
the corresponding circulation for various field configura- 
tions, and in particular for a quadrupole trap. In this 
connection we also note that the effects of a topological 
Aharonov-Casher phase in a toroidal configuration was 
studied earlier by Petrosyan and You || . 

Since very strong electric fields are required for the ef- 
fects of slac t° be non-negligible, it is important to inves- 
tigate how electric polarization effects could influence the 
potential trapping the atoms, and this is done in Sec. IV. 
Finally, we summarize our conclusions in Sec. V. 



II. EFFECTIVE VECTOR POTENTIALS IN THE 
GROSS-PITAEVSKII EQUATION 

We start with the Hamiltonian for a single particle with 
mass m, charge q, spin s and magnetic moment fi in an 
external electric and magnetic field, 



H 



1 

2m 



[p - gA(r) - a A c(r)} 2 - qA (r) - ft ■ B(r), (6) 



where — (A , A) is the usual EM gauge field, slac — 
fl x E/c 2 , and fl = //S. The spin operator S is a matrix 
acting on the (2s + l)-component wave function. For a 
spin 1/2 particle, and to linear order in the fields, this 
Hamiltonian follows from the Dirac equation (see, e.g., 
Q ) , and we will simply assume the same form of H for 
particles with general charge and spin, and in particular 
for neutral bosons. 

The second-quantized description of a collection of 
such fully-polarized neutral bosons with spin s is given by 
a coherent-state path integral for the 2s + 1-component 
complex scalar field ip with a Lagrangian, 



Ira 4 



(7) 



where we have added a contact-interaction term with 
strength A. 

The next step is to use an adiabatic approximation to 
"freeze out" the spin degree of freedom, assuming that 
the B field is strong enough to polarize all the atoms 
in a state given by the magnetic quantum number m s . 
Technically, we proceed by decomposing %j) — ipx(s) = 
?AU(s)xj, where if) is a single component complex scalar 
field, and x(s) a 2s + l component spinor satisfying x(s)^- 
x(s) = 1 and x(s)t n Sx(s) = s |10|. Here xt is a constant 
spinor corresponding to having the spin fully polarized 
in the, say, z-direction, and U is a unitary operator that 
rotates \] to the space time dependent spinor x(s(r,t)) 
that describes a spin at the position r pointing in the 
direction of the unit vector s(r). 

The adiabatic approximation for the spin dynamics 
now amounts to fixing this unit vector along the local di- 
rection of the magnetic field, i.e. taking s(r, t) — B(r, i), 
and neglecting fluctuations in the spin direction. (Note 
that we do not assume the orbital motion to be adia- 
batic.) The resulting path integral 

C = ^(ihd t + m s fJ,B)i/} - ^—ip\~ihV - a A c - &b) 2 ^ 

Am 



(8) 



is now over the single component field ip, which however 
couples to the effective gauge potential 



sl b = ih{U^VU) 



TT- 



(9) 
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It is easy to verify that the vector potential a#, is noth- 
ing but the Berry potential discussed in the introduc- 
ction. Also note that (for m s > 0) the term m s pB 
acts as a trapping potential for atoms. In the case of 
a time dependent B, there is also a time component, 
a*g, in the Berry potential, which would be relevant e.g. 
in a TOP trap. We again stress that the gauge poten- 
tials emerge because of the adiabatic assumption, which 
in this context amounts to ignoring spin-flip processes 
(which are suppressed by the Zeeman gap, assumed to 
be large enough.) 

The saddle point of the action of Eq. (|J) is determined 
by the Gross-Pietaevskii equation 



t;2 

— (- 



ihV - a eff ) 2 



m s pB - ^(<l>\ 



0, (10) 



■ is the condensate wave function, and a c ff = as - 



iTidt<j) '- 
where > 

We note that Eq. (|l^) resembles the time-dependent 
Ginzburg-Landau equation describing a superconductor, 
in the sense that there is a complex condensate wave 
function coupled to a space-dependent gauge potential. 
From this analogy we would expect tha an effective mag- 
netic field would either be expelled, or penetrate in the 
form of quantized vortices. We shall elaborate on this in 
the next section. 



III. EFFECTIVE POTENTIALS, CURRENTS 
AND CIRCULATION 

To understand the effects of the effective vector poten- 
tials &b and slaci we follow the corresponding analysis 
for a superconductor. There is, however, an important 
difference in that (at least in the approximation consid- 
ered here) the gauge field is not dynamical but rather a 
background field determined by the orbital motion. In 
particular there are no kinetic terms, and thus no reason 
to minimze the action with respect to the gauge poten- 
tials. ^From this follows that there can be no screening of 
the effective Berry or Aharonov-Casher fields, and con- 
sequently no Meissner effect. We thus expect our system 
to act as an extreme type II superconductor. (Only if 
kinetic terms together with h 2 B or h 2 AC terms were gen- 
erated by e.g. a renormalization group procedure, could 
there be a Meissner effect.) In this sense our system is 
also very similar to a bucket of helium II where rotation 
with an angular frequency u) can be described by a min- 
imal coupling to an effective gauge field h c g — ImQ. In 
our case, however, the effective field can have a more gen- 
eral space time dependence, just as in a superconductor. 

With these comments, we now proceed to minimize the 
ground state energy for a fixed effective magnetic field, 
and as in the cases of a superconductor or a rotating 
Bose-Einstein condensate, the essence of the argument is 
related to the single-valuedness of the phase of <f>. (For 
simplicity we consider only the zero temperature case. 



It is not hard to generalize to finite temperatures well 
below the Zeeman energy by introducing a suitable free 
energy.) 

For a homogeneous system, we parametrize <f> = 
^j~p§~e %s l h , with po being the constant mean atomic den- 
sity, and j = pqv = po(VS — a c s)/m the corresponding 
current density. The single-valuedness of the wave func- 
tion then gives the usual condition <f> dr- (mv + a G ff ) = nh 
on the circulation of the velocity field. For a cylindrically- 
symmetric geometry, and a circular path of radius r 
around the z axis, we get 



r d9[mve{z,r) + (a B + a A c)e] = nh. 



(11) 



As a first illustration, let us assume a c ff(r, 6, z) = 
o- c s(r)0 and take as an ansatz solution a vortex along 
the z axis with vorticity n. Equation (nil) then implies 



mrv(r) = nh- — $ s g(r), 

Z7T 



(12) 



where $ c ff(r) = $ac + $b is the flux due to the 
Aharonov-Casher and the Berry magnetic fields through 
the surface spanned by a circle of radius r. The kinetic 
energy per unit length e(n) associated with the vortex 
state is then given by 



(13) 



z(n) — J d 2 r—mp(r)v 2 
f R dr 

oc / — [nh- $ eff (r)] 2 , 
J£ r 



where £ is the healing length and p(r) is approximated 
by po for £ < r < R, and taken as zero elsewhere. 

Just as a rotating bucket of Helium, or a type II su- 
perconductor, we expect the presence of the magnetic 
flux to induce rotation in the form of vortex states. To 
get a rough estimate of when this happens, we assume 
a constant density and <E> e ffM = kr l . One can then 
show that for vortices to be energetically favourable, i.e., 
e(n) < e(0), one must have 



dr • a cff = $ e ff(i?) ~ nh, 



(14) 



where R is the radius of the condensate. Neglecting the 
Berry phase, this is nothing more than the condition 
that the Aharonov-Casher phase be of order 1, in agree- 
ment with the estimate of Ref. [0, where a toroidal 
geometry was used and a purely topological Aharonov- 
Casher phase was considered. The general case with non- 
constant density must be treated by solving the Gross- 
Pitaevskii equation, which can be done numerically, but 
we will not do it in the present study. 

Next we investigate the field configurations needed 
to achieve $ c ff ~ h. As a first illustrative (but com- 
pletely unrealistic) example, we consider a constant mag- 
netic field B = BqZ, and an electric field with a con- 
stant gradient, E = (E y/R)x, implying ag = and 
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a AC = (^E /R)y. Therefore, the effective magnetic 
field is constant, \)ac — (l^Eo/R)z. Thus we have 
$ e ff = irR 2 bAc = ttR^Eq/c 2 , which yields the condition 



eE$R ~ m e c 



(15) 



We now turn to a more realistic example, and assume 
that the atoms are trapped in a quadrupole field B = 
B'(x, y, —2z) by being in a state of polarization 



-/<-■ 



lie + yy — 2zz 



y 



-4z 2 



(16) 



i.e., opposite to the polarizing field B. This state is the 
so-called low-field seeking state (the absence of trapping 
at the center of the trap is a problem, since the atoms 
which are located there escape from the trap, and there 
are numerous tricks for fixing that, like in TOP traps or 
Ioffc-Pritchard traps; we do not worry about these extra 
complications in the present study.) 

With E = E z we have (at r — R, i.e., in the periphery 
of the cloud), 



1 . „ fiE R - 

a AC = ^xE= 9, 
c 2 c 2 VR 2 + 4z 2 



and 



$ac = 27I>£o 



R 2 



VR 2 + 4z 2 



(17) 



(18) 



The polarization also adds to a e ff the contribution from 
the Berry potential Eq. M): 



&b = hm s (cos (f — 1)9, 



hm s (cosip — 1) = —hm s 



2z 



VR 2 + 4z 2 



(19) 



(20) 



where cos <p = B z = —2zj\jR 2 + z 2 . This expression for 
the Berry potential is derived within a spin one-half rep- 
resentation by using the unitary operator U = e - i ^ m -< T / 2 ) 
where m = B x z/|B x z| describes the rotation axis, and 
ip = arccos(B • z) is the rotation angle. This potential 
has a Dirac string singularity on the positive z axis. In 
Fig. |l] we illustrate the corresponding fields. 

Since the particles are bosons, m s is an integer and the 
term —hm s in <&b corresponds to an integer number of 
flux quanta and induces no vorticity. The remaining part 
vanishes when integrated over the volume occupied by 
the entire cloud. Therefore, the Berry potential induces 
no vorticity. 

For z = 0, $4c (R) — 2-kR^Eq/c 2 , and we retain the 
result of Eq. ([l5]) obtained above. In fact, this result 
holds for general flux configurations of the form $ c g (r) cx 



/ / / 
I ' ' 



\ \ 



FIG. 1. The Aharonov-Casher field (higher) and the Berry 
field (lower) in polar coordinates. The vectors are normalized 
to be of constant length. 

The knowledge of the maximal electric field strength 
that can be applied to the condensate yields the minimum 
size of the cloud for our method of inducing rotation to 
be applicable. 



IV. LIMITS ON THE ELECTRIC FIELDS DUE 
TO THE ATOM POLARIZ ABILITY. 

Application of an electric field results in a Stark shift 
A, which in the quadratic regime can be thought of as 
the interaction between the electric field and an induced 
electric dipole moment, which is 



A = -d-E = --aE 2 . 

2 



(21) 



Here a is the polarizability of the atoms, and d = aE/2 
is the electric dipole moment. For alkali atoms a/fi ~ 
100 kHz/(kV/cm) 2 . The induced dipole moment changes 
the physics of the trapped atoms in several ways. The 
most direct effect is that the Stark shift changes the 
ground-state energy so that spatially-varying E fields 
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change the trap geometry. In addition, the dipole-dipole 
interaction given by 



Vdd(R) 



di-d 2 -3(di-R)(d 2 -R) 

AneR 3 



(22) 



where R is the distance between two atoms, and e is the 
permittivity, induces spin-flip transitions of the atoms 
into non-trapped states so they can escape. The dipole- 
dipole interactions do not conserve the spin but the mag- 
netic dipole-dipole interaction is sufficiently small (com- 
pared to the Coulomb interaction) for the spin relaxation 
rate to be low enough for the condensate to be observ- 
able. 

We have not attempted to calculate these quite com- 
plicated effects, but we made some estimates based on 
the following simple considerations. First, notice that 
the interaction given by Eq. ( p2[ ) differs from the mag- 
netic dipole-dipole interaction only by a constant fac- 
tor. In vacuum the interactions are of the same order 
when d/y/ea ~ A*\/M0' For these values the spin re- 
laxation rate should therefore be equally small. In the 
regime of the quadratic Stark effect, d = a.E/2, with 
a/h = 100kHz/(kV/cm) 2 , which implies E ~ 10 7 V/m. 

A more thorough study of the electric dipole-dipole 
interaction can be found in p"l) , where an effective scat- 
tering length is defined based on the effective potential 



V cS (R) = u Q 5 CR) + V dd (R) 



(23) 



where Uq = Anh a sc /m, with a sc being the scattering 
length for atom-atom elastic collisions described by the 
contact-interaction term in Eq. (|23]) . For 87 Rb, the ef- 
fective scattering length corresponding to V e g (R) is rel- 
atively unaffected up to E ~ 10 7 V/m. For higher values 
it changes drastically, attaining a negative value for 



£»6x 10 7 V/m. 



(24) 



This value actually provides an upper bound for E. We 
therefore estimate E ~ 10 7 V/m to be the upper limit 
for the field strength. Inserting this into Eq. (|l5|) yields 
that the lowest possible value of R is ~ 5 cm. 



V. SUMMARY 

We have demonstrated that inhomogeneous electric 
and magnetic fields applied in trapped atomic Bose- 
Einstein condensates can in principle be used in order to 
induce vortex states. More specifically, we have inferred 
the corresponding Gross-Pitaevskii equation satisfied by 
the condensate wave function in the presence of an elec- 
tric and a magnetic field, and we have investigated the 
influence of effective vector potentials such as the Berry 
and Aharonov-Casher type on the condensate wavefunc- 
tion. The Aharonov-Casher field can be used to induce 
vorticity in the cloud and we have found that the mag- 
nitude of the electric field E and the condensate size R 



have to satisfy the equation eER ~ m e c 2 — 0.5 MeV. 
As a relativistic effect, it requires large electric fields for 
typical sizes of the condensate. The magnitude of the 
electric field is in turn limited by the polarizability of the 
atoms. Estimating these limitations we find a lower limit 
for the condensate size to be R > 5 cm. 

Finally we want to stress that the mere possibility to 
study gauge interaction in Bose-Einstein condensates is 
of considerable interest. For instance, we know that con- 
densed systems in the presence of magnetic fields can 
exhibit different Meissner phases. As already indicated, 
this could be of relevance if the effective gauge fields 
somehow would aquire dynamics. Two dimensional fermi 
systems in strong magnetic fields exhibit various quan- 
tum Hall phases, and with this in mind it would be inter- 
esting to apply the methods developed in this paper both 
to e.g. lower dimensional atomic Bose-Einstein conden- 
sates and degenerate fermion systems. 
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